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Abstract 

It is shown that the Benney system is equivalent to a one body 
classical dynamical problem in which the interaction potential is a 
functional of the action function. The general solution of this equation 
is responsible for the general solution of the Benney system. Some par- 
ticular solutions of the selfconsistent Hamilton- Jacobi equation arising 
in this investigation are presented in explicit form. 







1 Introduction 



We use the form of the Benney system from [Q, where it is presented ('in 
its original form') as: 

dv ^ dv ( r (I dv ^dh ^ ^ ( ^ \ g (1) 

dt dx \Jo dx J dy dx ' dt dx \Jo J 

where the unknowns are the functions v = v{t; x,y), h = h{t, x). 

This paper is not aimed at the problems of symmetry and numerous rec- 
curence relations which follow from the Benney system (|I|) and are described 
in 0, H, ^, ^, 1^, H, but directly to the problem of construction of the general 
solution to it. The history and literature on Benney system reader can be 
found in 0. 

The general solution in our approach to the problem can provisionally be 
divided into two steps. In the first we propose a parametric representation of 
the solution of the Benney system (in implicit form) in terms of an assumed 
distribution function. In the second step we propose a self- consistent equa- 
tion, which this function should satisfy. The general solution of this equation 
is responsible for the general solution of the Benney system. 

Thus the goal of this paper is to demonstrate that the Benney system can 
be presented in the form of one body classical dynamic Hamilton- Jacobi equa- 
tion, with the potential expressed nonlocally in terms of the action function. 
We discuss also the connection of this problem with an ordinary differential 
equation of the second order. A particular solution of (|l|), depending on two 
arbitrary functions each of one argument is constructed. However the gen- 
eral solution of this equation must depend on one arbitrary function of two 
independent arguments and one function of one argument (the initial values 
of the functions v, h). Such a solution is unknown for us at this moment. 

The paper is organized in the following way. In section 2 we consider in 
detail the case with h = and explain the strategy of the calculations, which 
are also applicable to the general case without any considerable change. In 
sections 3 and 4 the general case is considered and shown that the Benney 
system is equivalent to a single non-linear integro-differentional equation, 
which in turn is equivalent to a Hamilton- Jacobi one, with the potential in 
the form of a nonlocal functional of the second order derivatives of the ac- 
tion function. In terms of the general solution (not only of its first integral!) 
of the self- consistent Hamilton- Jacobi equation the general solution of the 
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Benney system (|I]) is represented in implicit form. Two different particu- 
lar solutions of this equation with functional arbitrariness of two arbitrary 
single argument functions are also presented here. The connection of the 
self-consistent Hamilton- Jacobi equation with an ordinary differential equa- 
tion of the second order and the ways of obtaining its general solution are 
discussed in section 5. In section 6 we represent our results in the form of a 
theorem. Section 7 is devoted to a discussion of the results of the paper and 
discussion of possible perspectives for future investigation. In the Appendix 
the reader will find the connection of the solution of the main equation of 
the second subsection of section 3 with the Hamilton- Jacobi equation in a 
field of force depending only upon the time variable. 



2 Preliminary manipulations 

This section contains detailed calculations leading to a parametrical repre- 
sentation of the general solution of the Benney system in implicit form under 
the simplest assumption that /i = in (|l]). The reader for whom the result 
is more interesting than its derivation can begin reading directly the para- 



metrical representation of the solution (|29D , which is checked independently 
in the next few lines. 

Let us define u{t,x,y) = Jq dy'v{t,x,y') and resolve the second equation 
of the system (|l|) as 

h = H,it,x), u{t-x,H.^ = ~Ht (2) 
After such a substitution, the system (0) takes the form: 
Ut^y + UyU^y-u^Uyy + H^^ = Q, u [t] X , H ^) = ^ H f , u{t; x,0) = (3) 

To have an experience in working with the last system of equations let 
us at first restrict ourselves to the particular solution for which H = Const 
(really this is an inessential restriction for the parametrical representation of 
the solution as can be seen from the results of the section 3). 

Thus we have to solve the single equation for the unknown function u 
with the boundary condition: 

ut,y + UyUxy - UxUyy = 0, u{t; X, O) = (4) 
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Let us present the solution of in the form (j){u;t, x,y) =const, con- 
sidering (p as an unknown function with four arguments. Direct calculations 
of the first and second order derivatives using the rules of differentiation of 
implicit functions leads to the result: 
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The system of equations, which the functions a, (3, A satisfy (as a consequence 
of their definition only) takes the form: 



(7) 



In addition directly from the main equation for them, which can be 
written in two equivalent forms, follows: 



{0, Mtu = {(P, A}xy, {(p, a}yu = >^'^{(P, ^}yu 



(8) 



where {X, Y}a,b = XaYf, — XbYa. Considering a, (3 as functions depending 
upon four independent arguments (p,X,x,t, we obtain for the last equation 

of e: 



(A,-AA„)K-A2(^),) 



(9) 
and rewrite (J^) as: 

(At - a\u)Px + A = «x + {K - (10) 

(Aa; - pXu) = (Ay - \Xu)f3x, {Xy - AA„)aA = (Xt - aXu) 
Combining the first equation of (|10]) with the last two we conclude that: 

a,. = A, a = 0t, f3 = Qx (11) 
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and after substitution of(|33D into (|^) we arrive at a linear equation for a 
single function 0: 

- + Ae,,A = et,x (12) 

The most straightforward way to obtain its solution is by differentiation 
with respect to the argument A. This leads it to: 

d d 

The general solution of the last equation is obvious: 

e= d\'{\ - \')G{x + X't, A'; (p) + \K{x, t; 0) + L{x, t; 0) 



The initial equation (jT^) connects the functions K and L by a gradient de- 
pendence and, finally, solution of (p!2|), depending on two arbitrary functions 
G, F (each with three independent arguments) takes the form: 

= dX'{X- \')G{x + \'t,X';^) + Ft- (13) 

In what follows we will denote the "moments" of the function G by the 
corresponding upper index: 

G'{x, t; 0, A) = dX'iX'YGix + X't, A'; 0) 



In this notation a and p functions take the form: 

a = {Fu - Gl) - A(F,, - G°), (3 = (F,, - Gl) - A(F,.,, - G°.) 

Substituting a, (3, X via derivatives of the function (f) from (H) we come to 
the following system of equations of the first order: 

0. = {F.t - Gl)<j)u - (F.. - Gl)<f)y, cPt = {Fu - Gl)<f)u - {F^t - G'^^y (14) 

Up to now we have used only one equation from the system ([T0|) which 
together with the second equation (||) leads to (0). Two remaining equations 
(p!0D have as their direct corollary the following system: 



Xx — {Fxt — G\)Xu — {Fxx — G^^)Xy, Xt — {Ftt — G\)Xu — {Ftx — Gt)Xy (15) 
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The systems (JT^) and (|T5D are selfconsistent in the sense of equahty of the 
second mixed derivatives and our nearest goal now is to extract the conse- 
quences, which follow from them. 

Let us consider in ( p!4D the function (f> depending on two sets of four 
arguments (f) = 6{x,t,X;u) = ip{x,t, X;y). Keeping in mind (|15]) we can 
transform (|l^) to two forms: 

Ox = {Fxt — Gl)9u, 9t = {Ftt — Gl)du, 

^x = -{F^x - Gl)tPy, = -{F.,t - 
with the obvious general solution: 

^ = e = P[u + Ft-G\,\) = Q{y- Fx + G^, A) (16) 

where P and Q are arbitrary functions of their two arguments. Resolving 
( p!6D with respect to the first arguments of P and Q functions, we obtain a 
system of two functional equations, which the functions 4> and A satisfy: 

u + Ft-G^ = R{\(t>), y-Fx + G' = S{X,<f>) (17) 

Here R and S are arbitrary functions of two arguments (inverses of P, Q). 

We recall that F = F{x,t;(p), G^{x,t, X,(p), G^{x,t, X, ([)) and for this 
reason (0) is a system of functional equations which determine the functions 
A, (p implicitly. By the rules of differentiation of implicit functions it may be 
seen that the functions A and defined by (p!6D,(p!7D satisfy the systems (0) 
and (p!5|). 

To show this and to find the additional condition which connects the func- 
tions R, S (arbitrary up to now) let us calculate derivatives of the functions 
A, (p. For this the following abriviations will be used 

A = + — Ft^^, B = Rx + G\, C = Sff, — G^^ + F^^tp, D = Sx — G^ 

and the linear system of equations from which the corresponding derivatives 
are determined takes the form: 

Fu - G] = A<p, + BXt, -Fxt + = CcPt + DA*, 

Ftx -Gl = A(f)x + BXx, -Fxx + G° = G(f)x + DX^, 
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1 = A0„ + BXu, = C(/)u + DXu, 

= A(t)y + BXy, 1 = C4)y + D Xy . 

One can verify that the hnear systems ([T^) and ([T5|) are satisfied and for A 
we obtain: 

, 0, _ B _ Rx + Gl 
0„ D Sx- Gl 

From the last expression, keeping in mind the definitions of and G^ func- 
tions it immediately follows that: 

XSx = —R\ 

The last equation may be resolved in terms of only one arbitrary function T 
of two arguments A and (p: 

R = XTx~T, S = -T, 

To come back to the solution of the initial system and express w as a 
function of the arguments (x, t, y) it is necessary to put (f) =constant in all 
formulae above. After such a substitution the functions G and F reduce to 
functions of two indepentent arguments and T to only one, A. 

Substituting the last expressions into (pTf ) and recalling the boundary 
condition [H =constant) from we obtain: 

u + Ft~G^ = XTx-T, y-F, + G^ = -Tx, ,u(x,t,0) = (18) 

It is not difficult to check that the function T can be included in G by 
the of substitution: 

G{xi,X2) G{xi, X2) + T^^,^^{x2), 

which is equivalent to equating T in (|TBp to zero. 
Thus finally (|l^) takes the form: 

w + Fi - = 0, y - F, + = 0, u{x, t, 0) = (19) 

Two first equations define parametrically (via A = —Uy = —v) u as func- 
tion of its three arguments. The third one can be considered as a boundary 
condition, but really this is an equation determining two two dimensional 
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functions F{x,t), \{x,t,0) = v. Indeed, substituting u = y = into ( |I9| ) we 
obtain: 

Ft-G\x,t,u) = 0, -F^ + G'°(x,t,zy) =0 (20) 

Comparing the second mixed partial derivatives of the function F and keeping 
in mind the definition of the "moments" G^, G° we obtain finally: 

(z/i - uu^)G{x + ut, u) =0 (21) 

It is easy to check, that the solution of Monge equation: 

J^t - i^T^x = 

is implicitly defined by the equation 

G{x + vt,v) = g = constant 

Below we rewrite (^) including the function F in the lower limit of the 
integrals and emphasise the functional dependence of all functions involved: 

/•—v{t,x,y) /•—v(t,x,y) 

u{t,x,y)- dXXG{x + Xt,X) =0 y+ dXG{x + Xt, X) = 

Ju{t,x) Ju{t,x) 

(22) 

J^t - ^^x = 



The equations ([2^ ) define implicitly the general solution of the Benney 
system in the case h = 0, depending on one function G of two independent 
arguments (which is in connected with the initial value of the function v). 



3 The general case of the Benney equation 

Now we consider the general case of the Benney system (|l]). All steps of 
the calculations will be the same as the previous with the additional terms 
connected with H^x in (0). 
Equation (^ now reads: 

Dets (pu (puu (put + Dets (py (pyy (py^ + H^^(pl = (23) 

V (py (puy <Pyt J \<Pu <Puy <Pux J 
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The definitions of the functions a, /3, A do not change and the equations 
which follow from (10) are preserved in form. 



As a corollary of ( PB| ) the modified main equation (^) is: 

(A, - AA„)(aA - \\^)x) +H,. = (24) 

The equation for the function B takes the form: 

(A, - XXu) (6, - Ae,,A + e^.A) + H,. = (25) 

Further calculations word for word repeat the same story of the previous 
section. The crucial point is the understanding that an arbitrary function of 
four arguments 9 = 9(t, x, A, 0) may be presented in the form: 

e = / rfA'(A-A')G'(t,x,A';0) + F(t,x,0) + A$(t,x,0) (26) 

where G = 0a,a- Indeed (^) is one possible form a Taylor series for the 
function G with respect to one of its independent arguments. The form (^) 
doesn't have any connection with equation which the function G satisfy. 

All formulae between and (pO|) preserve their form with the obvious 
exchange —Fx—>-^,Ft^F and lead to the final result; instead of (^0]): 



^ + ^-^^ = 0, 7y + $ + G° = 0, u{x,t,0) = 0, u{x,t,H^) = -Ht, (27) 

where F, $ are two arbitrary functions of three arguments [t, x, <p). 

After differentiation of the second equation (^) with respect to operator 
4 — A-^ we obtain: 

ay au 

1 

{\y — Wu) - 



e 



Substituting the last expression into (psD and subsequent differentiation of 



the result with the respect to the argument A we obtain an equation for 
function G = 0a,a: 

-Gt + AG, - H,,Gx = (28) 



The last equation is exactly the equation for the distribution function for a 

1 

2' 



one body classical dynamical system with Hamiltonian function + H, 
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Indeed, the condition of conservation of some quantity along the trajectory 
can be written as: 

Pt = {H, p}a,x 

This is exactly the equation for the function G above. ( In the quantum 
case it is the equation for the density matrix with the exchange of Poisson 
brackets with commutators). 

Now let us include functions F, $ from ( pTj) into the first boundary con- 
dition and rewrite (|27|) in the form: 



u - [ d\XG{t,x,X) = y + [ dXG{t,x, X) = (29) 

Ju J u 

where is a function of two arguments t, xli is obvious that to satisfy the 
first boundary condition v = —v{t,x,0). [] 

The condition on the second boundary leads to the additional restiction: 

Ht + r dXXG{t, x,X) = r dXG{t, x,X) = Q (30) 

(yu = -v{x,t, -H^)). 

Now let us check by direct calculation the conditions under which func- 
tions u{t,x,y), H{t,x) defined by the ( pi?D and (pOj) are a solution of the 
Benney equation (||). 

In the notation used above the Benney equation may be rewritten as 
{uy = v = -A): 

Vt + vv^ - u^Vy + H^^ = 
After differentiation of the second equation (|29D with respect to t we find: 

P — V 

—VtG{t,x,—v) — utG{t,x,u) + / dXGt{t, X, X) = 



Substituting the derivatives of Gt from the equation for it (^Sp we obtain 
finally: 

/-"rfAAG,(i,i.A) - H„(G(-v) - Gil')) - i^tG)") 



^We remind the reader that in all formulae below it is necessary to put (j) ^constant, 
as explained in the previous section. 
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By the same technique we have: 

d\XG,{t, X, A) + {vv,G{-v) - vpM^)) 



UxVy 



G{-v) 



or 

J'" dXXG^{t, X, A) - yvxG{p) 

-{VV^ - UxVy) = — 

G{-v) 

After summation of these results we obtain the equahty: 

G{u 
G{^v) 



vt + vv^ - u^Vy = -H^x + — -{J^t - - Hx 



and conclude that in order to satisfy the Benney equation u function must 
be the solution of the equation: 

— lyUx — H^x = 

The condition of equality of the second mixed partial derivatives with 
respect to t, x of the function H follows from (|30D and leads to the conclusion 
that the function /i satisfies the same equation as i/ and the condition of 
selfconsistency of the whole construction my be written in the form: 

Hxx = ^t- T^T^x = fJ't- fJ-fJ-x, Hx= dXG{t, X, A) (31) 

Ju 

where G satisfies (^81) . 

Considering A in equation (^) as an unknown function of three indepen- 
dent arguments {x,t,g) and defining: 

g = G{x,t,X) 

we arrive at the equation, which this function satisfies: 

Xt - AA, -H,x = (32) 
from which it is follows that: 

ly = X{x,t,gi), fi = X{x,t,g2) 
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After exchaging integration variables under the integral sign in pO| ) we amal- 
gamate the Benney equation with the corresponding boundary conditions 
into a single (seemingly very strange) integro-differential equation: 



r92 

Xt - XXx- / gdgXg^^ = 



(33) 



The main function G from (p8|) may be reconstructed implicitly via the so- 
lution of the last equation: 



A = X{x,t,G) 



3.1 The simplest particular solution of the Benney sys- 
tem 

In this subsection we would like to demonstrate the selfconsistency of the 
previous construction by the example of the "simplest" solution of Benney 
system. This also will give us the possibility of understanding which class of 
functions is typical for this problem. 

Of course, distribution function G = —-^ = constant satisfies (|^) with 
an arbitrary choice of the potential function H^. Thus as a corollary of (^) 
we have {G = — 2z)- 

From (|29|) we obtain: 

v = -2Ay~u, u = -Ay^~iyy, h = = ^{fi - u) (34) 

By direct computation one can become convinced that (|3^) is really the 
solution of the Benney system if /i, u functions satisfy (|35|) , which in the case 
under consideration takes the form 

fJ-t - fJ-fJ-x = l^t - I'l'x = j^ifJ' - T^)x (35) 

For linearization of the last system let us exchange the dependent and 
independent variables. In other words let us consider (t, functions of 
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the pair of variables Corresponding formulae for this transformation 

are obvious: 



X = x{fi, v) = XyVt + x^nt, 1 = XyV^ + x^n^ 

After solving of the last equations and substitution of the result into (??) we 
obtain the first order linear system of equations: 

(x + vt)^ = -(x + fit)^ ^ ^'^^'^ + ^'"^ 
From the last system we conclude that 

{x + ut) = -6i,, {x + fit)=9^, 

t^'-^. .^-d^.he„-^e,), (36) 

fi — u — 1/ 2 K 

and obtain the equation, which the function 6 satisfy. We present it in terms 
of the variables S = |(// + z/), i? = — i/) in which it takes the most simple 
form ^ 

(1 - ^)6'e,s = Or,r 

This is a linear equation of second order with the separable variables. It 
allows us to obtain its general solution depending on two arbitrary functions, 
each of one argument. Equations ( p6| ) have to be reversed and variables /i, 
may be expressed as functions of the initial variables x, t. 

Formulae give a paricular solution of the Benney system depending 
on two arbitrary functions each of one argument. Of course, this is not the 
general solution of this system which must depend on one function of two 
arguments ( the initial value v{x,y,0)) and one function of one argument 
(the initial value h{x,0)). 

The example considered shows that it is impossible to expect some simple 
analytical expression for the general solution of the Benney system. 
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3.2 Particular solution of main equation for A function 

Let us seek the solution of main equation ( p3D in the form: 

X = xA{t,g) + B{t,g) 

We have in consequence: 

ixAt + Bt)-AixA + B) + i dgg{xAg + Bg)), = Q (37) 

Equating coefficients at zero and unity, decreasing x to zero and solving the 
system of equations arising for the functions A, B we obtain finally: 

A = _^±1R±XM ^^^^ ^gg^ ^ ^^g^^ (3g) 

t + U{g) Jg, 

where U{g)^V{g) are arbitrary functions. 

Resolving (implicitely) the equation (|38|) with respect to the function G 

A = \{x,t,G) 

and substituting the result into ( pOj) we come to the solution of the Benney 
system depending on two arbitrary functions f/, V , each of one argument. Of 
course this solution is of a different kind to that in the previous subsection. 
In this case the distribution function G is different from a constant. In 
some cases under a definite choice of the form U, V functions all calculations 
may be done explicitely and solution of Benney system can be presented by 
analytical expressions. 

In the Appendix the reader can find another way to obtain the solution 
of the present subsection, connected with the ordinary differential equation 
of the second order (see section 5). 

4 One body classical mechanical problem 
with self-consistent potential of interaction 

Equation (|33D has the form of two dimensional conservation law. It allows 
us to solve this equation in the form: 

A^ [92 

A = -'S'x, — + / gdgXg = -St 

Z Jgi 
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After eliminating the function A we pass to the one body Hamilton- Jacobi 
equation: 

St + ^- gdgSg,x = St + ^ + V{x, t) = 0, (39) 

where the potential function V is expressed nonlocally via the derivatives of 
action function S. 

We would like to point out here that a situation of this kind is not a new 
one. It is sufficient to remember the famous non-linear Schrodinger equation, 
where the potential of interaction coincides with the square of the modulus of 
the wave function. In the problem considered here the action function plays 
the same role (classical limit of the wave function). The potential function 
is expressed in a non-local way in terms of the derivatives of this (|39D. 

5 Connection with a second order ordinary 
differential equation 

After differentiation of the main equation (|33| ) with respect to the argument 
g (keeping in mind that the integral term of equation doesn't depend upon 
g), we obtain: 

(Ag)* = {\Xg)x, Ag = fx, XXg = ft 

where the function / satisfies in turn the equation arising after eliminating 
the function A from the previous system: 

i~p~)g ~ fx 

Jx 

Considering in the last equation x = X as an unknown function of the 
independent arguments (t, /, g) we have in consequence: 

f = f{X,t,g), l = fxXf, = ft + fxXt, = fg + fxXg (40) 

and as a corollary the equation, which the function X satisfies: 

X.Xtj - XfXt,, = 1 (41) 

Differentiation of the last equation with respect to the argument t leads to: 

^ttfXg — XttgXf = 
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This means that the function Xu depends on the arguments {g, f) only via 
the function X and thus we have 



Xu = Q{X,t) (42) 



In other words, from pT|) it follows that function X is the solution of a 
ordinary differential equation of the second order (|4^) with respect to the 
argument t. 

Let us present the general solution of this equation in the form: 

X = X{cuC2;t) 

and fix the choice of integration ci, C2 constants by the condition 0: 

Then C\ , C2 as functions of the arguments (/, g) thanks to pl|) are connected 
by a canonical transformation with the generating function W{ci^g): 

C2 = Wd{ci,g), f = Wg{ci,g) 



Collecting all these results we are able to rewrite the main equation (|33D in 
the form ofa functional integral equation determining the unknown function 
X: 

X = X{W,,,Ci, t), Xe.Xt,,, - X,,Xt,c, = 1, 

[92 

Xu{Wc„ci,t)-{ dggWg),, = 0, c2 = W,,{ci,g) (43) 

where, as was mentioned above, W{ci,g) is the generating function of the 
canonical transformation. 

In (^3]) the unknowns are both X and W and we have no idea at this 
moment of any constructive solution of this system, except by guesswork. In 
the Appendix we show that the choice of the second order ordinary differential 
equation in the form Xu = Q{t) corresponds to a solution of subsection 2 in 
the section 3. 



^Let us recall that these constants are defined up to a general covariant transformation. 
This circumstance always allows us to satisfy the condition below. 
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6 The main Theorem 

The results of the previous sections may be summarised in the following 
Theorem 

Each solution of nonlinear integro-differential equation 



r92 

Xt - AA^ - / dggXg^^ = 

-'51 



is connected (implicitly) with the solution of Benney system 



r dy—\ — + — - — + —(['' dyv] - 
io dx ) dy dx ' dt dx \Jo J 



dv ^ dv 

dt dx Wo dx J dy ' dx 



via the following formulae 

/•—v{y,x,t) /•—v(y,x,t) 

u{x,y,t)- dXXG{t,x,X) =0 y+ dXG{t,x,X)=0 

Ju{x,t) Ju(x,t) 

u{x,y,t)= dyv{y,x,t), h= dXG{t,x,X)= dggXg 
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where 

fi = X{t,x,g2), ly = X{t,x,gi), X = X{t,x,G) 



7 Outlook 

The main result of the paper consists in the Theorem of the previous section, 
giving asolution of the Benney system parametrically in implicit form in 
terms ofa distribution function and a selfconsistent equation (^31), which this 
function must satisfy. If it is possible to find the general solution of the last 
equation in the future the problem of the construction of the general solution 
of Benney system would be solved with the help of formulae (pOD. 

The most interesting questions about possible representations of the re- 
sults in the form available for experimental application (the origin of the 
Benney system is the hydro dynamical problem of surfaces waves [0]) are out 
of the framework of the present paper. To encompass these questions will 
be possible after more detailed investigation of the main equation (|33|) and 
the discovery of ways of its regular integration (maybe on the level of the 
computer computations). 
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All results of the present paper are obtained on intutive background and 
demand for the their rigorous foundation more powerfuU mathematical meth- 
ods. The authors can guess that they connected with the better known in- 
vestigation of the inner symmetry group of Benney system. We hope to come 
back to these interesting questions in our future publications. 
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9 Appendix 

In this Appendix we would try to widen the solution of the second subsection 
and obtain it on a more systematic background. With this purpose let us 
consider more precisely and evaluate the integral in (^). 




In the process of the evalution above we have used equations connected with 
derivatives of the functions / and and (0). In (0) the derivatives 

with respect to the argument t is only a partial one. To have the possibility 
of taking it out of the integration sign of it is necessary to increase it to a 
total derivative. This is achieved by the following obvious manipulation: 



Let us consider the possibility that the second term of the last equality in 
its turn is the total derivatives with respect to the argument t. It is obvious 
that for this it is sufficient assume additionally: 





(44) 





U{f,g) + a{t,g)f + b{t,g) 



(45) 



( but not the t variable). 
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Under the assumptions above we can once integrate (^) with the result: 
Xt + \ r dggait, g) = Pig, /), X = + P{g, f)t + Qig, f) (46) 

Z Jgi 

from equations ( PD and (|45|) we obtain additionally: 

Xf = {-{aJ + k))-"^, Xg = {-{aJ + bt))--^iU + af + b)) (47) 

The most simple way to resolve the last system of equations is to consider as 
a first step the condition of equality of the second mixed partial derivatives 
and compare one of the (selfconsistent) equations ( ^Tj ) with (^). 

The condition of selfconsistency of the second mixed derivatives reads as: 

atj + htg = at{U + af + b)- 2{U} + a){atf + h^) (48) 

In next three lines is presented the result of consequent differentiation of (^) 
with the respect to the argument /: 

atg = -at{Uf + a) - 2Uff{atf + h) 

= 3at-^ + 2{atf + bt) 
^fff 

= «*(3(^)/ + 2) 

Going in the reverse direction we obtain in consequence ( all big letters are 
arbitrary functions of argument g only): 

Uff = B{A-2f)-l, Uf = B{A-2f)-"2+D, U = -B{A-2fY^ + D f + E 

atg = -{D + a)at, 2b + Aa = F 

The last equation for the function a after integration once with respect to 
the argument t may go over to the form {a + D = z) of a Riccati equation 
(now with respect to the argument g: 

^. + Y = ^. + y 
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In the last equation we have represented the arbitrary function arising in 
a special form allowing us to integrate the Riccati equation with the final 
result: 

a + D = — — + —r—^ — 

After these preliminary calculations we are able to compare the first equation 
with (|46|) . For [atf + bt) from the results above we have: 

{aJ + h) = {A-2f)- 



or 

(Xf) -iPft + Qf) -(2737^$; 

All other computatrions are obvious and lead finally precisely to the solution 
of the second subsection in section 3. 
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